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We show that the size of the O(fl^) flavour violating cutoff artifacts that have been found to affect 
the value of the neutral pion mass in simulations with maximally twisted Wilson fermions is 
controlled by a continuum QCD quantity that is fairly large and is determined by the dynamical 
mechanism of spontaneous chiral symmetry breaking. One can argue that the neutral pion mass 
is the only physical quantity blurred by such cutoff effects. O(fl^) corrections of this kind are also 
present in standard Wilson fermion simulations, but they can either affect the determination of 
the pion mass or be shifted from the latter to other observables, depending on the way the critical 
mass is evaluated. 
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1. Introduction and main results 

Numerical data for the mass of the neutral pion in maximally twisted Lattice QCD (Mtm- 
LQCD) [|I|] simulations show unnaturally large lattice artifacts ' , despite the fact that on general 
ground they are expected to be 0{a^) corrections ||]. This is in striking contrast with the small- 
ness of the cutoff effects observed not only in the mass of the charged pions, which are related 
through the Goldstone theorem to exactly conserved lattice currents but also in all the other 
so far measured hadronic observables. Quite remarkably small lattice artifacts are found even in 
matrix elements where the neutral pion is involved [ ]T0| , [TT]]. 



In this talk, relying on arguments based on the Symanzik analysis [12] of lattice cutoff effects, 
we give an explanation for the origin of such peculiar corrections, showing that they are a general 
feature of any Wilson-Uke fermion regularization, whether twisted or not, and their appearance in 
the pion mass or instead in other observables depends on the choice of the twisted angle (zero or 
n/2) and the way the critical mass is determined. 

In sect. ^ we recall the properties of the Symanzik approach for the description of cutoff effects 
in LQCD with Wilson fermions and we discuss how the critical mass is determined. In sect. || we 
illustrate the nature of O(a^) artifacts in Mtm-LQCD and in sect. ^ how they show up in the standard 
Wilson fermion regularization. We end with some concluding remarks in sect. ^. 

2. Symanzik expansion and critical mass in Wilson fermion LQCD 

A) We consider Nf = 2 LQCD with quarks regularized as Wilson fermions. For generic values 
of the bare (twisted, /x, and untwisted, mo) mass parameters the lattice action reads 



~ a ^ la -> 

7- V - - V* V + csw jCJ • F + mo + /Airs-l^ 



X. (2.1) 



where for the sake of generality we have also introduced the clover term. In this talk we are 
interested in two specific cases comprised in ( |2. 1[ ). 



Mtm-LQCD, which is obtained from (p7l|), by setting ji = □(a") and mo = Ml^, where 



is some estimate of the critical mass. The physical interpretation of this regularization is most 
transparent in the so-called "physical basis", resulting from the field transformation 



i//- = exp(/7r75TV4)z, = ;i:exp(/7r75TV4) ^ (2.2) 

\lf. (2.3) 



7-V-/75t3 ( -^V*V + cswja-F+M,^,. ) +At 



• The clover standard Wilson fermions action, S'l, which is obtained by setting = and 
mo = m + M^j. with m an 0{a^) quantity. With this choice the most appropriate basis for discussing 



physics is the ;^ -basis itself in which eq. ( |2.1| ) was written in the first place 



'The neutral to charged pion mass splitting measured in the unquenched Mtm-LQCD simulations carried out in 
ref. with the tree-level improved Symanzik gauge action turns out to be smaller (and of opposite sign) than the 
quenched result ||^ where the standard plaquette gauge action was used. This finding is interesting in view of the 
established relation |Q ^ between the magnitude of this splitting and the strength of metastabilities detected in the 
theory at much too coarse lattice spacings kv. 
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B) The Symanzik effective Lagrangian associated to the Wilson LQCD action (2.1) reads 



^4 = ^™ + z[D + m + /75T>]Z : 



(2.4) 
(2.5) 



where the four-dimensional operator (all the necessary logarithmic factors are understood) specifies 
the continuum theory in which correlators are evaluated. The very definition of effective action, as 
a tool to describe the a dependence of lattice correlators, implies that the mass parameters in ^4, if 
not exactly vanishing, must be 0(a°) quantities. Thus all the lattice artifacts affecting M^^. will be 
described by operators of the form a'^A'^j^xX^ ^ = 1,2, ... in 5^sym- 

After using the equations of motion of ^4, the 0(a) piece of 5^sym reads (^5;sw and 5i are 
0(1) coefficients) 

^5 = ^5;SwZ'f^ • Fx + 5iA^cD^;i: + 0{m,n) . (2.6) 

The terms multiplied by powers of m and/or /i are not specified in eq. ( p^ ) because they are not of 
relevance for the topic discussed in this note. We recall that the coefficient ^5;sw vanishes, if csw 



in eq. (2.1) is set to the value appropriate for Symanzik 0(a) improvement. 
The O(a^) part of d^sym has a more comphcated expression of the type 

3 14 
^6 = I%<r+^6;4^7M(£»M)';^ + I^6;,-<J>6;,- + 52A^CD^Z + 0(m,At), (2.7) 



(=1 



i=5 



where the first three operators are purely gluonic, the fourth is a chiral (but not Lorentz) invariant 
fermionic bilinear and the remaining ones are four fermion operators, which we find useful to write 
in the form (equivalence with the list in can be proved using Fierz rearrangement) 



^6,5 = 


izzMx), 


<J>6;6 = 


<J>6;7 = 


-{XY5X){XY5X) , 


<J>6;8 = 


<J>6;9 = 


{XYkX){XYkX) , 


<J>6;10 = 


<I>6;11 = 


-- {XYxY5X)iXYkY5X) , 


<J>6;12 = 


<J>6;13 = 


-- {XOxvX){X(yxvX) , 


<J>6;14 = 



i:b{x^'x){x^'x), 

-Lh{XY5r''x)iXY5^''x), 
'Lh{XYi^^X){XYk^''x), 



(2.8) 



C) Both in the case of standard Wilson and twisted mass fermions the condition to determine 
the critical mass is the vanishing of the PCAC mass. Let us examine these two cases separately. 
1) tm-LQCD - The condition (it is convenient to rotate the quark fields to the i//-basis (|2.2D) 



«^L((v^)'otV)(^,0(W5tV)(o)) 







(2.9) 



leads to a determination of the critical mass which is "optimal" {M°f^) in the sense that with this 
choice all the leading chirally enhanced cutoff effects are eliminated from lattice correlators |^. 



In the spirit of the Symanzik approach the condition ( |2.9| ) must be viewed as a relation holding 
true parametrically for generic values of a (and jj.). As a result, it is equivalent to an infinite 
set of equations, each equation corresponding to the vanishing of the coefficient of the term a'^, 
k = 0,l,2,.... >From the vanishing of the a° term one gets 



d^x { ( YYo y) {x,t){ VY5 T V) (0) ) 



0, 



(2.10) 
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by which restoration of parity and isospin is enforced. This means that, if (the 0{a^) piece of) 



mo is chosen so as to verify eq. ( |2.9| ), then we will simultaneously have m = in ( |2.5| ) and the 
identification on the lattice of i/ZyoT^VA with the (time component of the) vector current (the 
identification of i/ZysT^i//^ with the pseudoscalar density being trivial). The further implications 
of eq. (23) are conveniently exposed looking at its Symanzik expansion. In ref. [^] it was proved 
that at 0(a) implies the condition 



cont 

^^'"^ = bs-sw a • F VA + 5i A^cD ViYs t V + 0(At ) • 



(2.11) 
(2.12) 



Eq. ([2.11[) should be read as a constraint fixing 5i. At 0{a ) the only relevant term (tni] is the 



one where VqP^ is inserted with (the integrated density) ^g^"". The latter in the i/A-basis has the 
expression 

' + ^A3jcDW5TV + 0(Ai^), (2.13) 



^^P— even 



with parity-even. Since in the continuum limit (because of parity invariance) one gets 

Sd^xSd^y{^^-''''^''{y)V^{x)P\id))\^ont = 0, the condition implied by (gj) yields % = 0, owing 
to / (Px J d^y{\j/i'Y5Z^\lfiy)Vo{x)P^{0))\cont 7^ 0. It follows from this analysis is that the estimate 
of the critical mass provided by ( ^ ) is not affected by O(a^) effects. These arguments can be 
generalized to all orders in a and show that MS'' can only display 0{a^'^^^), = 0, 1, ... corrections. 
The latter are determined by constraints, like ( 2.11 ), that fix the value of the coefficients d2p+i in 
front of Yijs '^^ W- 

2) Standard clover Wilson fermions - The condition for the vanishing of the PCAC mass is 



3oD(A[;(x,o/'''(o)) 



2LAP''i^,t)p^m 



■ 'MPCAC 



@n=0, 



(2.14) 



which apart from the normalization and a trivial time derivative is exactly eq. ( [2.9| ) (though written 
in the ;^ -basis) with the only difference that now jj. =0. This condition is in practice implemented 
by looking for the limiting value of mo for which mpcAC 0^- 

The vanishing of the twisted mass is at the origin of all the differences resulting from the two 
ways of subtracting the Wilson term. In fact, if /i is set to zero, from the symmetries of the Wilson 
theory and the associated Symanzik expansion one cannot conclude anymore that 82 vanishes. 



Rather at 0{a ) eq. ( |2.14| ) fixes the value of ^2 through the condition 



{KiO)\^f\KiO)) 



0, 



(2.15) 



where is the full six-dimensional operator of the Symanzik Lagrangian associated to the clover 
improved Wilson fermion regularization (including the contribution of the two fermion operator 
a52AQ^Y)XX)- Iri general discretization errors of any order in a will affect the critical mass deter- 



mination (2.14) (except those Unear in a owing to clover improvement). 
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3. Neutral and charged pion mass in Mtm-LQCD 



• Neutral pion mass - The quantity of interest for the study of the neutral pion mass is the 
zero-momentum four-dimensional Fourier transform of the two-point (subtracted) correlator 



n{p)=a'Y^e'P^{P\x)P\Q)) 



(3.1) 



It is immediate to recognize that at /? = and in the limit of very small lattice pion mass one gets 

|2 



rL(o) 



(3.2) 



>From the Symanzik expansion of Ti{0) through orders a^ included one can prove [^] that, even 
in the absence of the clover term, thanks to the optimal choice of the critical mass (see sect. §), one 
arrives at the equation 



\G.r 



mt 



1 



mt 



+ 0(- 



mt 



(3.3) 



where the continuum pion mass has been simply indicated by mjt. Consistently with the results of 
;^PT 15, 4[ [l^, |l^], a simple Taylor expansion leads to the key formulae of this note 

ml + a^t;, + 0{a\x\ y) , l^, ^ {k' (0) lif^" \^\^))\c 



I cont 1 



^Mtm 



<^P— even 
■^6 



. (3.4) 



• Estimating O(a^) lattice artifacts in \i - To estimate the size of the O(a^) artifacts (| 
we need to compute l^ji in the chiral limit. This can be done under the assumption that a sufficiently 
accurate estimate of l^j^, can be obtained in the vacuum saturation approximation (VSA). Quenched 
studies show that VSA works quite well for matrix elements of four-fermion operators between 



pseudo-scalar states OlSQ. We must then identify the operators in ( |2.8D that have non-vanishing ma- 
trix elements between states as nij^ — > and give a non-zero contribution in the VSA. An exam- 
ple is P^P^ = (vZ/sT-'i/a) (i/z/jT-^i/a) which corresponds in the hst ( |Z8| ) to the operator —{xx) iXX)- 
Noticeably one can prove that the matrix elements between states of the four-fermion operators 
in ^g'^'™ of interest for our evaluation of ^j,: are all proportional to (7r^(0)[/'^P^|7r^(0))|cont in the 
limit m| 0. Thus up to a numerical factor in the VSA we can write a^l^j^ a'^\Gji\'^, with Gji the 



continuum (renormalized) analog of the quantity defined in ( |3.2[ ). 

An estimate of Gj^ can be obtained either by a direct lattice measurement of a^Gji iQ] or 
exploiting the WTI 2mq{Q.\P^\n^)\com = fn^i^- Using the results of |Q |Tl]], the two evaluations 
turn out to be numerically well consistent yielding jG^^P ~ (570 MeV)"^, a number ~ 20 - 25 times 
larger than the typical scale Aq(-.q ~ (250 MeV)'*. 

• Charged pion mass - Replacing the isospin index 3 in eq. ( |3.l| ) with either 1 or 2, one finds 

:m2±+a2(;r±(0)|^^^'-|;r±(0))|eont + O(aV^,a4) 



■■ m\ + 0(a^m?.,fl:'*) 



(3.5) 



The last equality follows from the invariance of under SU(2)ob = {Q\,Q\,Qv) and it is in 

perfect agreement with ;^PT [Q, 16, 17] and the similar result derived in ref. [Q. The lattice square 
pion mass splitting in Mtm-LQCD can thus be estimated up to terms of 0{a^m\,a'^) with the result 



Am, 



Mtm 






L 




L 



■a^C^r ~a^(570 MeV)'' 



(140MeV)2 @ a^^ ~2.3GeV. 



(3.6) 
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This number compares very nicely with the value of the splitting (180(40) MeV)^ reported in iQ]. 

• Where else does l^ji enter? - Given the impact we have seen it has on the lattice expression 
of the neutral pion mass, an important question to ask is where else (besides m^j \i and all related 
energy factors) can the key parameter 1^^ appear in the Symanzik expansion of lattice quantities. 



The answer requires a detailed analysis which we have no space to report here Jl IQ. The outcome 
of it is that to all practical purposes the only interesting place where l^n enters is just m^j \i. 

4. Pion mass and O(a^) artifacts with standard Wilson fermions 

Proceeding as before, one gets for the (clover improved) standard Wilson pions the formula 

ml\L = ml + a\%[0)\^i\%{0))l,,, + O{aW^,a^), (4.1) 

where isospin indexes are understood owing to the SU(2) flavour symmetry of the lattice Wilson 
theory. No 0(a) terms are present {b^-sw = 5i = in (p^)) as we are assuming clover improvement. 

The particular way in which the critical mass is fixed reflects itself into the form of the 
Symanzik effective Lagrangian of the lattice theory. Here we will examine two choices. The 



first corresponds to the standard procedure where the critical mass is determined from eq. ( |2. 14| ). 
The second is somewhat more exotic and corresponds to fixing the critical mass by using the deter- 
mination provided by Mtm-LQCD. 



• The standard way affixing the critical mass - At O(a^) the condition for the vanishing of 



'MpcacIl implies the relation (2.15), which fixes &i in terms of other parameters of the theory and 



in particular of the matrix elements of the four-fermion operators (2.8) between one-pion states 



Given the uniqueness of the Symanzik expansion (eqs. (2.4), (2.5)) for Wilson fermions (close 
to the chiral limit), a theoretical analysis similar to that we have sketched in sect. ||, together with 
the numerical estimate of G^, shows that there are sizable contributions in eq. ( |2.15| ). As a result 
chances are that ^2 ^ 1 because this coefficient has to compensate for the large value of l^ji. The 
consequences of this situation are twofold. 

1) No O(a^) artifacts will affect the value of the lattice pion mass because they are absent in 
'MpcAC \l thanks to ( 2.15 ) and, as it follows by taking the limit f ^ oo in eq. ( 2.14 ), one has 

2 fn 



I -niYCAC ■ (4.2) 

/ICj^tl L L 

2) On the contrary, in other observables, like for instance the mass of the hadron h, there will 
appear O(a^) effects proportional to {h\xx\h) multiphed by the possibly large coefficient ^2. 

• Using the critical mass of Mtm-LQCD - If the estimate of the critical mass as determined 
in Mtm-LQCD is instead employed, since, as shown in [^, only odd powers of a come into play, 
the term o^I^qq^xX will not appear in the Symanzik Lagrangian. As a result one will have ^2=0 



and the 0(a ) corrections to the square pion mass in eq. (4.1) will not be zero. But now, no 0(a ) 



corrections stemming from a^h^QQY)XX will affect other observables. 

5. Concluding remarks 

In this talk we have argued that there are peculiar O(a^) cutoff effects in LQCD with Wilson 
fermions which have a dynamical origin related to the mechanism of spontaneous chiral symmetry 
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breaking. In Mtm-LQCD they only affect the neutral pion mass making it substantially different 

from that of the charged pion. If the standard Wilson fermion regularization is employed, where 
these discretization errors will show up will depend on the way the critical mass is determined. With 
the usual determination, pion masses are free from these lattice artifacts, but the latter will appear 
in other physical quantities, such as hadronic masses. If, instead, the critical mass as determined in 
tm-LQCD is employed such O(a^) terms will only affect the value of the pion mass. 
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